Relativistic Green's function approach to charged-current neutrino-nucleus 

quasielastic scattering 



Andrea Meucci, Carlotta Giusti, and Franco Davide Pacati 



CO 

o 

O 

CN 

> 

o 

CN 



oo 
o 



cn 

o 



o 

=5 



X 



Dipartimento di Fisica Nucleare e Teorica, Universita degli Studi di Pavia and 
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, F27100 Pavia, Italy 
(Dated: February 8, 2008) 

A relativistic Green's function approach to inclusive quasielastic charged-current neutrino-nucleus 
scattering is developed. The components of the hadron tensor are written in terms of the single- 
particle Green's function, which is expanded on the eigenfunctions of the nuclear optical potential, so 
that final state interactions are accounted for by means of a complex optical potential but without 
a loss of flux. Results for the (v^fjT) reaction on ie O and 12 C target nuclei are presented and 
discussed. A reasonable agreement of the flux-averaged cross section on 12 C with experimental data 
is achieved. 

PACS numbers: 25.30.Pt: Neutrino scattering, 13.15.+g: Neutrino interactions, 24.10.Jv: Relativistic mod- 
els, 24.10. Cn: Many-body theory 



I. INTRODUCTION 

The reactions with an incident neutrino interacting 
with a complex nucleus have gained in recent years a 
wide interest, owing both to astrophysical reasons and 
to the aim of investigating the neutrino properties with 
a high accuracy. Besides the measurements with large 
underground detectors, some experiments have also been 
performed (jj- L3- L3l using a pion beam which weakly de- 
cays producing leptons. In this case the most part of the 
neutrinos which are obtained is related to muons, with a 
smaller component of electron neutrinos. 

Both weak neutral- and charged-current scattering 
have stimulated detailed investigations. In particular, 
we are here interested in charged-current reactions at an 
energy below 1 GeV as they have shown to be depen- 
dent on nuclear structure effects. Different approaches 
have been applied to investigate such processes, includ- 
ing the so-called "elementary particle model" Q , random 
phase approximation (RPA) in the framework of a rela- 
tivistic Fermi gas model || or a Fermi gas model with 
local density approximation @, shell model 0, and 
relativistic shell model El, R PA among quasiparticles Q 
and continuum RPA |lfTlllj | . The reaction goes through 
a quasielastic (QE) mechanism, where the neutrino inter- 
acts with one single neutron and a proton together with 
a negative muon are emitted. The effect of final state 
interactions (FSI) has been stressed to significantly con- 
tribute to the cross section an d has been calculated 
with a relativistic shell model including a phenomcnolog- 
ical optical potential, which describes the interaction of 
the outgoing nucleon with the residual nucleus, with and 
without imaginary part [l3l | . 

The optical potential is fitted to reproduce the elas- 
tic proton- nucleus scattering through its real component, 
while its imaginary part takes into account the scatter- 
ing towards the inelastic channels. This means that the 
reaction channels are globally described by a loss of flux 
produced by the imaginary part of the complex poten- 



tial. This model has been applied with great success to 
exclusive QE electron scattering ^4|, where it is able to 
explain the experimental cross sections of one-nucleon 
knockout reactions in a range of nuclei from carbon to 
lead. In an inclusive process, however, where some of 
the reaction products are not detected and the inelastic 
channels are also included in the experimental cross sec- 
tion, the flux must be conserved. This fact is sometimes 
described by dropping the imaginary part of the optical 
potential. This procedure conserves the flux but it is not 
consistent with the exclusive reaction, which can only be 
reproduced with a careful treatment of the optical poten- 
tial, including both real and imaginary parts |l4j . 

In this paper we apply a Green's function approach 
where the conservation of flux is preserved and FSI are 
treated in the inclusive reaction consistently with the ex- 
clusive one. This method was discussed in a nonrelativis- 
tic and in a relativistic ^(| framework for the case 
of inclusive electron scattering and it is here adapted, 
in a relativistic framework, to charged-current neutrino 
scattering. In this approach, the components of the nu- 
clear response are written in terms of the single-particle 
optical-model Green's function. This result was origi- 
nally derived by arguments based on the multiple scatter- 
ing theory [l^ an( A successively by means of the Feshbach 
projection operator formalism 0, 0, 0, |2(j • The spec- 
tral representation of the single-particle Green's function, 
based on a biorthogonal expansion in terms of the eigen- 
functions of the non-Hcrmitian optical potential, allows 
one to perform explicit calculations and to treat FSI con- 
sistently in the inclusive and in the exclusive reactions. 
Important and peculiar effects are given in the inclusive 
(e, e') reaction by the imaginary part of the optical po- 
tential, which is responsible for the redistribution of the 
strength among different channels. 

In Sec. II the general formalism of the charged-current 
neutrino-nucleus scattering is given. In Sec. Ill, the 
Green's function approach is briefly reviewed. In Sec. 
IV, the results obtained on 16 O and 12 C target nuclei are 
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presented and discussed. Some conclusions are drawn in 
Sec. V. 



II. THE INCLUSIVE CROSS SECTION 

In a charged-current process neutrinos and antineutri- 
nos interact with nuclei via the exchange of weak- vector 
bosons and charged leptons are produced in the final 
state. The cross section of an inclusive reaction where an 
incident neutrino or antineutrino, with four-momentum 
/cf = (ei,ki), is absorbed by a nucleus and only the out- 
going lepton, with four-momentum /c M = (e, k) and mass 
to, is detected, is given by the contraction between the 
lepton tensor and the hadron tensor, i.e., 

(1) 



da = 2, IT W» v 



2 (2tt) 3 ' 

where G ~ 1.16639 x 1CP 11 MeV -2 is the Fermi constant 
and i? c is the Cabibbo angle (cos$ c ~ 0.9749). 
The lepton tensor is 

L „u = un^l T 7 5 )^ ± 7 5 )7^/, (2) 

spin 

where the upper (lower) sign corresponds to neutrino (an- 
tineutrino) scattering. After projecting into the initial 
neutrino (antineutrino) and the final lepton state, one 
has 

-^Tr [ 7 • k r (1 T 7 5 ) 7 • h t"] , (3) 

Z£j£ 



which can be written, by separating the symmetrical and 
antisymmetrical components, as 

1 " " (4) 



Sit 



l^TlT], 



with 



fcj* k v + k» fe M 



where ^ va!i 

^0123 _ 



gT h ■ k 
■0, (5) 
is the antisymmetric tensor with €0123 — 

e""~ = 1. 

Assuming the reference frame where the z-axis is par- 
allel to the momentum transfer q — ki — k and the y- 
axis is parallel to ki x k, the symmetrical components 
I a 1% V Vg ', and the antisymmetrical ones I a', I a 1 , I j^, as 
well as those obtained from them by exchanging their 
indices, vanish. 

The hadron tensor is given by bilinear products of the 
transition matrix elements of the nuclear weak charged- 
current operator J M between the initial state | ^o) of the 
nucleus, of energy Eq, and the final states | \?f), of energy 
Ef, both eigenstates of the (A+ l)-body Hamiltonian H, 
as 



^(q) I *o> 



*o | | * f ) S(E 



Et). 



(6) 



and involves an average over initial states and a sum over 
the undetected final states. The sum runs over the scat- 
tering states corresponding to all of the allowed asymp- 
totic configurations and includes possible discrete states. 

The transition matrix elements are calculated in the 
first order perturbation theory and in the impulse ap- 
proximation, i.e., the incident neutrino interacts with 
only one nucleon while the other ones behave as specta- 
tors. The current operator is assumed to be adequately 
described as the sum of single-nucleon currents, corre- 
sponding to the weak charged current 



f = [ F Y(Q 2 )r+i—FY(Q 2 )a^q v 
- G a (Q 2 )7' 1 7 5 + ^p(Q 2 )^7 5 ]t ± , 



(7) 



where r are the isospin operators, k is the anoma- 
lous part of the magnetic moment, q^ = (u>,q), with 
uj 2 , is the four-momentum transfer, and 



and F2 are the isovector Dirac 



Q = \q\ 2 
a* v = (i/2) [7^, 7" 
and Pauli nucleoli form factors, which are taken from 
Ref. [2l|. Ga and Fp are the axial and induced pseu- 
doscalar form factors, which are usually parametrized as 



G A 



9A 



(1 
2MG A 
ml + Q 



Q 2 /M\f 



(8) 
(9) 



where <?a = 1-267, is the pion mass, and Ma — 1.032 
GeV is the axial mass. 

The most general covariant form of the hadron tensor is 
obtained in terms of the two independent four- vectors of 
the problem, i.e. the four-momentum transfer q^ and the 
four-momentum P M of the target. The symmetrical and 
antisymmetrical components of the tensor can therefore 
be written as 



W» v 



w% 



Wi + W 2 q^q u + W 3 P^P V 
Wa (PV + q»P") 

W 5 (P»q v - q»P v ) + W 6 e^q a P f3 , (10) 



where the invariant form factors Wi are functions of the 
two scalars which can be formed with q^ and P M , i.e., 
Q 2 and P ■ q. From Eq. illOfl . it is clear that in our 
reference frame W 0x , W 0y , W xz , W yz , and W x s v vanish 
together with the tensor components obtained from them 
by exchanging their indices. 

The inclusive cross section for the QE z/(i/)-nucleus 
scattering, obtained from the contraction between the 
lepton and hadron tensors, can therefore be written as 



de df2 47r 2 



wo-Roo + v zz R zz 
vozRoz + v T R T ± v xy R xy . (11) 
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The coefficients v, obtained from the lepton tensor com- 
ponents, are 

vq = 1 + k cos d , 
v zz = 1 + fccosz? - 2 ] l , ' sin 2 # , 



voz = y- [1 + kcostf) + 



2 

q\e 



= 1 — K COS U H : — rr— sin V 



y^jz = — j — |— (l — fccos^) — -j— j— 
9 ; \q\e 



(12) 



where fc = \k\/e, $ is the lepton scattering angle, and 
m the mass of the emitted lepton, e.g., the muon mass 
~ 105.9 MeV. It was taken 



i* = ^(ir+<F). 



(13) 



taking advantage of the fact that W xx = W vv , as can be 
deduced from Eq. I|10(l . 

The response functions are given in terms of the com- 
ponents of the hadron tensor as 



Roo 


= w oa , 




Rzz 


= w zz , 




Rqz 


= w 0z + 


w z0 , 


Rt 


_ w xx + 


W yy ; 


Rxy 


= i(W yx 


_ w xy ) 



(14) 



III. THE RELATIVISTIC GREEN'S FUNCTION 
APPROACH 

All nuclear structure information is contained in the 
response functions and therefore in the hadron tensor. 
We apply here to the inclusive QE y{u)-mxc\eus scatter- 
ing the same relativistic approach [l|| which was rather 
successful in reproducing the cross sections of the in- 
clusive QE electron scattering. Here we recall only the 
main features of the model. More details can be found 
in Refs. [3 El 13 

For the inclusive process the components of the hadron 
tensor in Eq. © can equivalently be expressed as 

W^(u,q) = -ilm<*o I rHq)G(Et)J^q) | *„) ,(15) 

in terms of the Green's function G(Ef) related to the 
nuclear Hamiltonian H, i.e., 



G(E { ) = 



1 



E { - H + it] 



(16) 



Here and in all the equations involving G the limit for 
r/ — > +0 is understood. It must be performed after calcu- 
lating the matrix elements between normalizable states. 



The current operator is assumed to be adequately de- 
scribed as the sum of the single-nucleon currents of Eq. 

It was shown in Refs. 0, 0, EE E3 for the inclu- 
sive (e, e') scattering that the components of the nuclear 
response in Eq. i|15|) can be written in terms of the single- 
particle Green's function G{E) 1 whose self-energy is the 
Feshbach's optical potential. An explicit calculation of 
the Green's function can be avoided by its spectral rep- 
resentation, which is based on a biorthogonal expansion 
in terms of the eigenfunctions of the non-Hermitian op- 
tical potential V, and of its Hermitian conjugate V^, i.e., 



[£-T-V\E)] \&\E)) = 0, 
[£-T-V{E)\\&\E)) = 0. 



(17) 



Note that E and £ are not necessarily the same. The 
spectral representation is 



/>OG 

G(E) = / d£\x { i\E)) 



1 



E-£ +ir) 



(18) 



The hadron tensor can be reduced to a single-particle 
expression and can be written in an expanded form as 



W^(uj,q) = --Vim / 

x rr(£,E f -e 



d£ E~^ 

m £jf — e n — t 



irj 
(19) 



where n denotes the eigenstate of the residual Hamil- 
tonian of A interacting nucleons related to the discrete 
eigenvalue e n and 

T^{£,E) = X n { Vn \f\q)^l-V{E)\x { s\E)) 

x ( X <f\E) | y/l-V(E)j«(q) | <p n ) . (20) 

A„ is the spectral strength of the hole state | ip n ), which is 
the normalized overlap integral between | ^o) and | n), 
while the factor y/T— V'(E) was shown in Ref. E3l to 
account for interference effects between different channels 
and to allow the replacement of the mean field V by the 
phcnomcnological optical potential Vl- 

After calculating the limit for r/ — > +0, Eq. I|19fl reads 



d£ 



M 



Ef — e„ — £ 



ReT% v (E { -s n ,E { -s n ) 



lmT^{£,E f -e n ) 



(21) 



where V denotes the principal value of the integral. Eq. 
(|21|l separately involves the real and imaginary parts of 
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The second matrix element in Eq. I(20() , with the inclu- 
sion of V%i an d disregarding the square root correction 
due to interference effects, is the transition amplitude for 
the single-nucleon knockout from a nucleus in the state 
| ^o) leaving the residual nucleus in the state | n). The 
attenuation of its strength, mathematically due to the 
imaginary part of V , is related to the flux lost towards 
the channels different from n. In the inclusive response 
this attenuation must be compensated by a correspond- 
ing gain, due to the flux lost, towards the channel n, by 
the other final states asymptotically originated by the 
channels different from n. In the description provided 
by the spectral representation of Eq. 1(181) . the compen- 
sation is performed by the first matrix element in the 
right hand side of Eq. 1)20(1 . where the imaginary part 
of V has the effect of increasing the strength. Similar 
considerations can be made, on the purely mathematical 
ground, for the integral of Eq. 1(21(1 . where the ampli- 
tudes involved in T£ v have no evident physical meaning 
when £ ^ E{ — e n . We want to stress here that in the 
Green's function approach it is just the imaginary part 
of V which accounts for the redistribution of the strength 
among different channels. 

The cross sections and the response functions of the 
inclusive QE neutrino (antineutrino)-nucleus scattering 
are calculated from the single-particle expression of the 
hadron tensor in Eq. ((21|) . After the replacement of the 
mean field V(E) by the empirical optical model potential 
Vl(-E), the matrix elements of the nuclear current oper- 
ator in Eq. I(2U|) . which represent the main ingredients 
of the calculation, are of the same kind as those giv- 
ing the transition amplitudes of the electron induced nu- 
cleon knockout reaction in the relativistic distorted wave 
impulse approximation (RDWIA) [24|]. Thus, the same 
treatment can be used which was successf ully app lied to 
describe exclusive (e,e'p) and ("f,p) data [HiHg. Here, 
of course, the nuclear electromagnetic current must be 
replaced by the nuclear weak charged-current operator 
of Eq. . 

The relativistic final wave function is written, as in 
Refs. m terms of its upper component following 

the direct Pauli reduction scheme, i.e., 

X { f\E)=[ 1 f+ * j , (22) 

V M + £ + S\E)-V^{E) P f+ / 

where S(E) and V(E) are the scalar and vector energy- 
dependent components of the relativistic optical po- 
tential for a nucleon with energy E [2fil ]. The upper 
component, 5'f+, is related to a two-component spinor, 
$t, which solves a Schrodinger-like equation contain- 
ing equivalent central and spin-orbit potentials, obtained 
from the relativistic scalar and vector potentials [53, , 
i.e., 

* f+ = \jD\{E)*t , (23) 

D £ (E) = 1+ SW-VW 
iV ' M + £ K ' 



where Dg(E) is the Darwin factor. 

As no relativistic optical potentials are available for 
the bound states, the wave functions ip n are taken as 
the Dirac-Hartree solutions of a relativistic Lagrangian 
containing scalar and vector potentials |29|, |30| . 



IV. RESULTS 

The calculations have been performed with the same 
bound state wave functions and optical potentials as in 

Refs. mum, 

where the RDWIA was able to repro- 
duce (e,e'p), (7,p), and (e, e') data. 

The relativistic bound state wave functions have been 
obtained from Ref. [2|j, where relativistic Hartree- 
Bogoliubov equations are solved in the context of a rel- 
ativistic mean field theory and reproduce single-particle 
properties of several spherical and deformed nuclei [3(| • 
The scattering state is calculated by means of the 
energy-dependent and A-dependcnt EDAD1 complex 
phenomenological optical potential of Ref. [26| , which is 
fitted to proton elastic scattering data on several nuclei 
in an energy range up to 1040 MeV. 

The initial states | tp n ) are taken as single-particle one- 
hole states in the target. A pure shell model is assumed 
for the nuclear structure, i.e., we take a unitary spectral 
strength for each single-particle state and the sum runs 
over all the occupied states. 

The results presented in the following contain the con- 
tribution of only the first term in Eq. 1(21(1 . The cal- 
culation of the second term, which requires integration 
over all the cigenfunctions of the continuum spectrum of 
the optical potential, is a complicate task. Its contribu- 
tion has been estimated to be small in the kinematics 
explored; hence, it is neglected in the present calcula- 
tions. 

In order to show up the effect of the optical poten- 
tial on the inclusive reaction, the results obtained in the 
present approach are compared with those given by dif- 
ferent approximations. 

In the simplest approach the optical potential is ne- 
glected, i.e., V = V 1 * = in Eq. JT7J|, and the plane 
wave approximation is assumed for the final state wave 
functions x^ - ' an d X ■ I n this plane wave impulse ap- 
proximation (PWIA) FSI between the outgoing nucleon 
and the residual nucleus are completely neglected. 

In another approach the imaginary part of the optical 
potential is neglected and only the real part is included. 
This approximation, that was sometimes used in the past, 
conserves the flux, but it is inconsistent with the exclu- 
sive process, where a complex optical potential must be 
used. Moreover, the use of a real optical potential is 
unsatisfactory from a theoretical point of view, since the 
optical potential has to be complex owing to the presence 
of open channels. 

The partial contribution given to the inclusive process 
by the sum of all the integrated exclusive reactions with 
one-nucleon emission is also shown in the following for a 
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FIG. 1: The cross sections of the 16 0(z^,,^i _ ) reaction, in- 
tegrated over the muon angle, for E v = 300, 500, and 1000 
MeV. Solid lines represent the result of the Green's function 
approach, dotted lines give PWIA, long-dashed lines show the 
result with a real optical potential, and dot-dashed lines the 
contribution of the integrated exclusive reactions with one- 
nucleon emission. Short dashed lines give the cross sections 
of the 16 0(z/ M , reaction calculated with the Green's func- 
tion approach. 



comparison. In this case only the eigenfunctions X °f 
V' are included and the imaginary part of the potential 
produces an absorption which does not conserve the total 
flux. We note that in the Green's function approach of 
Eqs. (|20[1 and i|21|) FSI are treated, consistently with the 
exclusive process, by means of a complex optical poten- 
tial. The imaginary part, however, does not produce a 
global loss of flux and it is responsible for the redistribu- 
tion of the strength among different channels. 

As a first study case we have considered the 16 O tar- 
get nucleus, for which the adopted single-particle wave 
functions and optical potentials have given a good agree- 
ment between RDWIA calculations and (e, e'p), (j,p), 
and (e, e') data. 

In Fig. \I\the cross sections of the 16 0(f M , reaction, 
integrated over the muon scattering angle, are displayed 
as a function of the muon kinetic energy T M for three dif- 
ferent values of the incident neutrino energy E v = 300, 




^ [deg] 



FIG. 2: The cross sections of the 6 0(v fi , fi~) reaction, inte- 
grated over the muon energy, for E v — 300, 500, and 1000 
MeV as a function of the scattering angle of the outgoing 
muon Line convention as in Fig. 



500 and 1000 MeV. The behaviour of the calculated cross 
sections is similar for the different energies. The effect of 
the optical potential increases with and decreases in- 
creasing E v . At 300 MeV, the result of the PWIA is 
much higher than the one of the Green's function ap- 
proach, while at 1 GeV the two results are almost the 
same but at the highest values of T M . The sum of the 
exclusive one-nucleon emission cross sections is always 
much smaller than the complete result. The difference 
indicates the relevance of inelastic channels and is due 
to the loss of flux produced by the absorptive imaginary 
part of the optical potential. In contrast, the cross sec- 
tions calculated with only the real part of the optical 
potential are practically the same as the ones obtained 
with the Green's function approach. Although the use 
of a complex optical potential is conceptually important 
from a theoretical point of view, the negligible differences 
given by the two results mean that the conservation of 
flux, that is fulfilled in both calculations, is the most im- 
portant condition in the present situation. In contrast, 
significant differences are obtained with a real optical po- 
tential in the inclusive electron scattering 0, 0] . 



6 



7 



- 
- 










f . , — • — ' * 


* / 
* / 


— " 




/' 


• / / 




•' / / 




■*/ 7 




:// 




/• s 




Is 






100 200 300 400 500 600 700 800 900 1000 

E„ [MeV] 



FIG. 3: The total cross section of the 16 0(i/ M ,/i~) reaction, 
integrated over the muon energy and angle, in terms of the 
neutrino energy E v . Line convention as in Fig. Q 

Qualitatively similar results are obtained in Fig. [21 
where the cross sections integrated over the muon energy 
are displayed for E v = 300, 500, and 1000 MeV as a 
function of the scattering angle of the outgoing muon. 

The global effect of FSI is clearly shown in Fig. [3] 
where the cross sections are integrated over the energy 
and the angle of the outgoing muon. The PWIA result 
is always much larger, while the loss of flux produced 
by the absorptive optical potential in the exclusive pro- 
cesses produces a too small cross section. In contrast, an 
optical potential with only a real component seems able 
to give a result comparable with the one of the Green's 
function approach. Small differences are found only at 
higher neutrino energies. 

In Figs. 2] |2 and [3] also the cross sections for the 
16 0(^ M , ) reaction are shown for a comparison. They 
are always much smaller than the corresponding cross 
sections with an incident neutrino. 

The different approaches have been compared with the 
experimental results of the LSND collaboration at Los 
Alamos for the 12 C(z^i, reaction y. L| . The cal- 
culations have been flux-averaged over the Los Alamos 
neutrino spectrum. In Fig. 21 the Green's function ap- 
proach is normalized to the experimental data, while the 
other results are accordingly scaled. The shape of experi- 
mental data is reasonably reproduced. The flux-averaged 



1 500 

1250- * 



1000- 




T M [MeV] 

FIG. 4: The distribution of the muon kinetic energy for the in- 
clusive 12 C(^ M , fi~) reaction. Experimental data from Ref. [J- 
The result of the Green's function approach is normalized to 
the experimental data. The other curves are scaled, accord- 
ingly. Line convention as in Fig. Q 

cross section integrated over the muon energy gives 11.15 
x 10~ 40 cm 2 . At the low energy of the experiment, how- 
ever, other processes, which are here not included, be- 
sides the quasielastic scattering can affect the inclusive 
reaction, in particular the excitation of the discrete states 
of 12 C. The experimental value (10.6 ± 0.3 ± 1.8) xlO~ 40 
cm 2 Q is slightly overestimated. The results obtained by 
other calculations 0,0,0 IS El El El gi ye larger values. 



V. SUMMARY AND CONCLUSIONS 

We have applied to and (^ A ,,// + ) reactions 

an approach based on the spectralization of the single- 
particle Green's function in terms of the eigenfunctions 
of the complex optical potential and of its Hermitian con- 
jugate. This approach has proved to be rather successful 
in describing inclusive electron scattering. Its advantage 
stands in the fact that it is able to include in a simple 
way and keeping flux conservation the final state inter- 
actions by using an optical potential which is essential to 
reproduce exclusive electron knockout reactions. 

The method is applied within a relativistic framework 
to weak charged-current reactions for an energy up to 1 
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GeV, where nuclear structure effects are important and 
dominant with respect to nucleon-resonance excitations. 
The reaction mechanism is assumed to be a direct one, 
where the incident neutrino interacts with only one neu- 
tron and a proton is emitted together with a negative- 
charged muon. A single-particle model is used to de- 
scribe the structure of the nucleus and a sum over all 
single-particle occupied states is performed. 

Calculations for the 16 O target nucleus have been pre- 
sented at neutrino energies up to 1 GeV. The optical 
potential and flux conservation have a large effect on the 
cross sections. For 12 C, the results averaged over the ex- 



perimental flux of neutrinos are compared with the avail- 
able data. A fair agreement is obtained in arbitrary units. 
The integrated cross section results somewhat larger than 
the experimental value. 



Acknowledgments 

We would like to thank W.M. Alberico for useful dis- 
cussions. 



[1] M. Albert, est al., Phys. Rev. C 51, 1065 (1995). 
[2] C. Athanassopoulos, et al, Phys. Rev. C 56, 2806 (1997). 
[3] L.B. Auerbach, et al, Phys. Rev. C 66, 015501 (2002). 
[4] S.L. Mintz and M. Pourkaviani, Nucl. Phys. A594, 346 
(1995). 

[5] H. Kim, J. Piekarewicz, and C.J. Horowitz, Phys. Rev. 

C 51, 2739 (1995). 
[6] S.K. Singh, N.C. Mukhopadhyay, and E. Oset, Phys. 

Rev. C 57, 2687 (1998). 
[7] A.C. Hayes and I.S. Towner, Phys. Rev. C 61, 044603 

(2000) . 

[8] C. Volpe, N. Auerbach, G. Colo, T. Suzuki, and N. Van 

Giai, Phys. Rev. C 62, 015501 (2000). 
[9] W.M. Alberico, M.B. Barbara, S.M. Bilenky, J.A. Ca- 

ballero, C. Giunti, C. Maieron, E. Moya de Guerra, and 

J.M. Udi'as, Nucl. Phys. A623, 471 (1997). 
[10] E. Kolbe, K. Langanke, and P. Vogel, Nucl. Phys. A652, 

91 (1999). 

[11] N. Jachowicz, K. Heyde, J. Ryckebusch, and S. Rom- 
bouts, Phys. Rev. C 65, 025501 (2002). 

[12] C. Bleve, G. Co', I. De Mitri, P. Bernardini, G. Mancar- 
ella, D. Martello, and A. Surdo, Astropart. Phys. 16, 145 

(2001) ; G. Co', C. Bleve, I. De Mitri, and D. Martello, 
Nucl. Phys. Proc. Suppl. 112, 210 (2002). 

[13] C. Maieron, M.C. Martinez, J.A. Caballero, and J.M. 
Udfas, Phys. Rev. C 68, 048501 (2003). 

[14] S. Boffi, C. Giusti, F. D. Pacati, and M. Radici, Elec- 
tromagnetic Response of Atomic Nuclei, Oxford Studies 
in Nuclear Physics, Vol. 20 (Clarendon Press, Oxford, 
1996); S. Boffi, C. Giusti, and F. D. Pacati, Phys. Rep. 
226, 1 (1993). 

[15] F. Capuzzi, C. Giusti, and F.D. Pacati, Nucl. Phys. 
A524, 681 (1991). 



[16] A. Meucci, F. Capuzzi, C. Giusti, and F.D. Pacati, Phys. 

Rev. C 67, 054601 (2003). 
[17] Y. Horikawa, F. Lenz, and N.C. Mukhopadhyay, Phys. 

Rev. C 22, 1680 (1980). 
[18] C.R. Chinn, A. Picklesimer, and J.W. Van Orden, Phys. 

Rev. C 40, 790 (1989); 40, 1159 (1989). 
[19] P.M. Boucher and J.W. Van Orden, Phys. Rev. C 43, 

582 (1991). 

[20] F. Capuzzi and C. Mahaux, Ann. Phys. (N.Y.) 254, 130 
(1997). 

[21] M.J. Musolf and T.W. Donnelly, Nucl. Phys. A546, 509 
(1992). 

[22] J.D. Walecka, in Muon Physics, Vol. II, edited by V.H. 
Hughes and C.S. Wu (Academic Press, New York, 1975), 
p. 113. 

[23] S. Boffi and F. Capuzzi, Nucl. Phys. A351, 219 (1981). 
[24] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. C 64, 
014604 (2001). 

[25] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. C 64, 
064615 (2001). 

[26] E.D. Cooper, S. Hama, B.C. Clark, and R.L. Mercer, 
Phys. Rev. C 47, 297 (1993). 

[27] B.C. Clark, in Proceedings of the Workshop on Rela- 
tivists Dynamics and Quark-Nuclear Physics, edited by 
M.B. Johnson and A. Picklesimer (John Wiley & Sons, 
New York, 1986), p. 302. 

[28] M. Hedayati-Poor, J.I. Johansson, and H.S. Sherif, Nucl. 
Phys. A593, 377 (1995). 

[29] W. Poschl, D. Vretenar, and P. Ring, Comput. Phys. 
Commun. 103, 217 (1997). 

[30] GA. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C 55, 
540 (1997). 



